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Abstract - We have developed a novel method to evaluate the potential profile of a molecular 
motor at each chemical state from only the probe's trajectory and applied it to a rotary molecular 
motor Fi-ATPase. By using this method, we could also obtain the information regarding the 
mechanochemical coupling and energetics. We demonstrate that the position-dependent transition 
of the chemical states is the key feature for the highly efficient free-energy transduction by Fi- 
ATPase. 



A molecular motor is a chemical engine that converts 
^^hcmical free energy into mechanical motion; to under- 
s tand the design principles of the molecular motor, it 
I— is essential to reveal its energetics. During the conver- 

sion, chemical potential is converted into mechanical po- 

^tential that then powers the motor's mechanical motion. 
\^ Thus far, for an ATP-driven rotary molecular motor, Fi- 
OOATPase, or Fi-motor [TJj^, the outputs of the energy 
transduction such as the heat dissipation through rota- 
^^tiions IlHS]j and its maximum work per step against ex- 
C\l ternal torque [7] have been measured. These studies have 
^~^r,evealed that the Fi-motor can transduce the free energy 
change of ATP hydrolysis to mechanical motion at nearly 
. . 100% efficiency. However, the mechanism to achieve such 
^ a high efficiency remains unclear. To reveal this, infor- 
^> mation about the mechanical potential that drives the 7- 
5^ shaft's rotation is crucial. 

5^ , In this letter, we describe our new method for recover- 
ing the mechanical potential profile of each chemical state 
only from the single-molecule trajectory (Fig. and its 
application to the rotational trajectory of the Fi-motor. 
By using this method, we could also obtain information 
regarding the mechanochemical coupling and energetics. 
The method is simple and applicable to not only other 
molecular motors but also systems in a broad range. 

Fi -motor is a rotary molecular motor, whose central 7 
shaft rotates uni-directionally against a stator ring (as/Ja) 
while hydrolyzing ATP to ADP and phosphate (P;) (Fig. 
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Fig. 1: Fi-motor (03/337 sub complex). The 7 shaft rotates 
unidirectionally when ATP is hydrolyzed to ADP and Pi. By 
attaching the 0:3/33 ring to a glass cover slip and a dimeric probe 
particle to the 7 shaft, the rotation can be observed under an 
optical microscope. 



[T^). At low ATP concentrations, Fi-motor exhibits a step- 
wise rotation with pauses corresponding to ATP waiting 
states. During this state, each /3 subunit has one of the 
following three different chemical states: with ATP, with 
ADP (and Pj), and empty [5111] (Fig. HJi). The chemical 
states of the three /3 subunits rotate cooperatively due to 
the ATP binding, ATP hydrolysis, and the release of ADP 
and a phosphate. Since the conformation of the /3 subunit 
varies with its chemical state [2l[T0], the interaction po- 
tential between the (3 subunits and 7 shaft rotates and 
drives the 7-shaft's 120° rotation. This picture leads to a 
simple model that the stepwise rotation is described as a 
Brownian motion on a potential corresponding to the ATP 
waiting state, wherein the potential's position shifts 120° 
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Fig. 2: Recovery of potentials from single- molecule trajectories, 
a, Simplified reaction scheme of Fi-motor. ot subunit is omitted 
for simplicity. The binding of ATP to the empty /3 subunit 
triggers 120° rotation, b, Interaction potentials between the 7 
shaft and /3 subunits. c, Single-molecule trajectory. We recover 
the potential profiles (b) from single- molecule trajectories (c). 



Table 1: List of symbols. 



L 


Number of video frames 


N 


Number of chemical states 


X 


Trajectory of the motor's position, 
(xi, . . . 


s 


Chemical-state sequence, (si, . . . , sl) 


s* 


argmaxCT(5, X). 

s 


a(X S) 


Score function to be maximized 




Path probability that X realizes under S 


a{S) 


Path probability that S realizes 




Probability that x realizes at s 


P{s) 


Probability that s realizes 


Ws{x'\x) 


Transition probability from x to x' at s 


W{s'\s) 


Transition probability from s to s' 


D 


Diffusion coefficient 


Us{x) 


Potential profile at s 


T 


Sampling period 



upon the ATP binding (Fig. [2Jd). Such discrete model is 
known to well describe the Fi-motor's stepwise rotation 

Analytical methods. — Let X and 5* denote col- 
lectively the observed trajectory of the motor's probe 
(a;i,...,XL) and the chemical-state sequence of the mo- 
tor (si, . . . , sl-i; Sk G {1, • ■ • , N}), respectively, where L 
is the number of the video frames and N is the number 
of states. In the case of Fi-motor at low ATP concentra- 
tions, = 3. The essence of our method is to find the 
most probable chemical-state sequence S* that maximizes 
a well-defined score function a{S,X) for the given X: 



S* = argmax<7(S', A) 

s 



(1) 



We define the score function as the joint path probability 
of A and S: a{X,S) = asiX)a{S). Here, 



L-2 



a(5) = P(si)[]W^(sk-hi|sk) 



(2) 



fe=i 



is the Markovian path probability that 5* realizes, P{s) is 
the probability that a state s realizes, and W^(s'|s) is the 
transition probability from a state s to another state s' . 
On the other hand, 



L-l 



CTs(A) = _Psi(a;i) Y\_ Ws^ix]^+i\xk) 



(3) 



k=l 



is the Markovian path probability that A realizes under a 
given S; here, Ps{x) is the probability of observing a; at a 
state s, and Ws{x'\x) is the transition probability from x 
to x' in successive two frames at s. Let D and r denote the 
diffusion coefficient and the sampling period, respectively. 
When the mean square displacement between successive 



frames Dt is sufficiently small, Ws{x'\x) becomes approx- 
imately 



1 



exp 



- [x' - X + ^d^Usix)^ /4Dt 



(4) 



to the first order of Dt [T2], where /cbT is the thermal 
energy and Us{x) is the mechanical potential at a state s. 

In the below simulations and experiments, we calculated 
D as the mean power spectrum density of the rotational 
rate in the frequency range from 0.8 kHz to 1.2 kHz since 
the fluctuation-dissipation relation is expected to hold at 
the high frequency region [5]. For Ps{x), we used Gaussian 
functions fitted to the peaks of the probe's angular distri- 
bution. For evaluating a(X,S), we need to know Us{x), 
P{s), and W(s'|s) additionally, which are not known a 
priori. We iteratively evaluate these unknown parameters 
as well as S* as follows. 

The algorithm consists of three steps. On the ba- 
sis of rough estimations of the parameters (Step 1), we 

find S* = argmax(T(X, 5*) by using Viterbi algorithm |13j 

s 

(Step 2). Then, we update Usix), P(s), and W{s'\s) on 
the basis of AT and S* (Step 3). We iterate Step 2 and Step 
3 alternatively typically 20 times. cr(A, S*) does not al- 
ways increase monotonically over iterations. We adopt S* , 
Us{x), P{s), and W^(s'|s) that maximize a(A, 5**) as the 
final estimation. The detailed procedures are as follows. 

Step 1: Initialization. As the initial estimation, we use 
U^{x) ^ -lnps(x), P{s) = 1/N, and W{s'\s) = M/L, 
where M is the total number of steps estimated from the 
trajectory. 

Step 2: State identification. We find S* that rigorously 
maximizes (t(A, S). This seems hopeless since the number 
of possible 5 is A^-^ (A^ = 3 and L ~ 2, 500, 000 in the 
experiment below). However, the Viterbi algorithm finds 
iS** at a computational cost of 0{N'^L), which depends on 
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L linearly |13] . 

Step 3: Potential estimation and parameter update. We 
divide the spatial coordinate to bins and calculate the 
number of transitions between the bins for every consec- 
utive two frames. Let be the number of transitions 
from the i-th bin to the j-th bin at a state. The equilib- 
rium distribution at this state, {tTi}, satisfies the balance 
condition, ttj = Y.i^i^i'^U where a;i_>j = rii^-J Y^k^'^^^- 
By solving this equation, we obtain the potential profile in 
this state as Ui = —k^TlwKi with a trivial additive con- 
stant. We fit Ui by using a sixth-order polynomial func- 
tion for estimating Us{x). We update P(s) and W^(s'|s) 
to P{s) = m^liL - 1) and W{s'\s) = rris^s' /iTT-s, where 
rus and mg-^s' are the numbers of s and steps from s to s' 
included in S**, respectively. 

Numerical simulation. — The validity of this 
method was tested by employing a simulated trajectory of 
a Brownian particle. The particle was assumed to be con- 
strained in a non-harmonic potential and obey a Langevin 
equation with parameters consistent with the below exper- 
iment (Fig. 13^). The potential's position was shifted with 
120° spacings randomly in a Poissonian process. Figure 
[3)3 demonstrates that our method recovered the poten- 
tial profile of each state accurately only from the stepwise 
trajectory. By using this method, we can also evaluate 
the angular positions where the transitions occur (6'tr); if 

7^ ^t+i ^ certain k, Xi^+i is the transition position 
for this transition. The distribution of 0tr was similar to 
the overall angular distribution (Fig. This correctly 
reflects the fact that we shifted the states randomly. The 
estimation error of 6'tr was 0.3 ±9.8° (mean ± SD). D was 
estimated to be 13.9 rad^/s in good agreement with the 
simulation parameter (13.8 rad^/s). 

Experiment. — We applied the method to the rota- 
tional trajectory of the Fi-motor at a low ATP concentra- 
tion (0.5 fj.M ATP, 2.5 /iM ADP, and 1 mM Pi) wherein 
the ATP binding was the rate-limiting step. The 7-shaft's 
rotation was probed by a 300 nm dimeric particle and ob- 
served at 9,000 Hz with an exposure time of 40 fis by using 
a high-speed camera (Fig. 2^) . Such short exposure times 
provide information regarding instantaneous angular po- 
sitions. Since backward steps are rarely observed in the 
absence of external loads [7] , we restricted the analysis to 
only the forward steps (ATP-hydrolytic direction). Note 
that our method can be applied to systems with bidirec- 
tional steps as well. 

Figure [IJd shows the potential profiles corresponding to 
the three ATP-waiting states recovered from the single- 
molecule trajectories. The potential had a simple struc- 
ture with a single minimum. The potentials are plotted 
with vertical shifts with the amount of the chemical free 
energy change of ATP hydrolysis (A/i ~ 16.5kBT, see Ma- 
terials and methods) . We found that the transitions do not 
occur at random positions but only at a limited angular 
region ahead of the potential minimum in the rotational 
direction (Fig. HJ:). It is noteworthy that the potentials 




Fig. 3: Numerical simulation of the Langevin equation, "fx — 

-d.u{x{t) - xo{t)) + m-, im) = 0; (mat)) = ^^k^Tm- 

Here, 7 = 0.3pNnms/rad^ is the rotational frictional coeffi- 
cient (D = 13.8rad^/s), T = 300 K, the frame rate is 9,000 
Hz, and the frame length is 2 x lO''. U{x) = O.Ss*^ - 0.5a:^ - 
3a::'* -f 3a::"^ -f lOx^ — 2x, where x is expressed in radians. The 
position of the potential xq is shifted with 120° spacings uni- 
directionally randomly in a Poissonian process at a rate of 6 
Hz. These parameters are chosen to be consistent with the 
experiment. Transitions within 20 ms after the previous tran- 
sition are prohibited since such successive transitions are not 
expected in experiments, a, Numerically generated trajectory 
coloured according to the estimated states, b, Recovered po- 
tential profiles (symbols) and the theoretical profile U (x) (solid 
lines), c. Angular distribution at transitions (histogram) and 
the overall angular distribution (solid line). The overall an- 
gular distribution is scaled by eye to be compared with the 
transition angular distributions. The potential profile in b and 
angular distribution at transitions from this state in c are de- 
picted using the same colour. 

intersect around this region. See Figs. [S^ and b for the dis- 
tributions of the mean transition angle ((^tr)) relative to 
the potential minimum (^btm) and the intersection point 
of the potentials (0x), respectively. 

It is expected that ATP frequently attacks the /3 subunit 
irrespective of the 7-shaft's angular position. However, 
our results suggest that the 7-shaft has to be ahead of the 
potential minimum for the ATP binding to proceed (Fig. 
[5j:); otherwise, the attack fails. Such a position-dependent 
substrate binding supports the "diffusion catch" mecha- 
nism for the torque generation [T^KH], in which the ther- 
mally activated conformation suitable for the substrate 
binding is stabilized by the binding. On the other hand, 
after the transition, the 7-shaft rotates along the poten- 
tial of the new state into its minimum driven by the bend- 
ing motion of the /3-subunit (power stroke). Thus, the 
mcchanochemical picture of the torque generation was elu- 
cidated in terms of the basic concepts only from the trajec- 
tory. Although we have assumed W to be independent of 
the position for simplicity, the estimated transition prob- 
ability depended on the position. It is straightforward to 
use a position-dependent W; this would increase the reli- 
ability of the estimation. 
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Fig. 4: Experiment. Recovered free-energy potentials and tlie 
angular distribution at transitions from the trajectory of an 
Fi-motor molecule, a, Rotational trajectory of the Fi-motor's 
probe particle observed at 9,000 Hz with an exposure time of 40 
/is. The trajectory comprises approximately 650,000 frames, b, 
Recovered potentials corresponding to the ATP-waiting states. 
The vertical lines denote the intersection points of the poten- 
tials. The potentials are plotted with vertical shifts with the 
amount of the chemical free energy change of ATP hydrolysis 
(A/i = 16.5 /ubT). c. Angular distribution at the transitions 
(histogram) and the overall angular distributions (solid line). 
See the caption of Fig. [3] for details. 

Furthermore, by using this method, we can evaluate 
thermodynamic quantities such as the heat dissipation of 
the probe only from the trajectory. The heat dissipated to 
the environment through the rotational degree of freedom 
of the probe per step, Q, is equal to the potential change 
between successive transitions [16]. Our result, (6'tr) — 6'x, 
implies 

Q^^^JL (5) 

(Fig. [S};). In fact, the estimated value of Q was 17.8 ± 
l.Qk-QT (mean ± SD) and was similar to A/x (16.5 A;bT). 
The slight difference between Q and A/i might be due to 
the estimation error of {On)- Because of the steep slope 
in the left-hand side of the potential (Fig. |4}d), a small 
estimation error of (0tr) results in a large estimation error 
of Q. It is possible that (^tr) estimated from the probe's 
motion is slightly different from that of the 7-shaft and bi- 
ased toward the potential minimum due to the time delay 
of the probe's motion. If this is the case, (0tr) is under- 
estimated and therefore Q is overestimated. On the other 
hand, because of this steep potential, when A/i is varied 
by controlling the concentrations of ATP, ADP, and Pi, 
the changes in 9^ are small. Therefore, the relation (O is 
expected to hold over a broad range of A/i provided rea- 
sonably that (0tr) does not depend on A/i largely. This 
suggests that Fi-motor transduces most A/i to rotational 
motions at a nearly 100% efficiency [S]. This is consis- 
tent with a previous study [S], wherein we had measured 
Q from the violation of the fluctuation-dissipation rela- 



tion using Harada-Sasa equality [T7|. It required the re- 
sponse measurement against an external load, which is 
accessible only in limited systems. The present method 
extracted a thermodynamic relation ([5]) only from the tra- 
jectory and furthermore demonstrated that the position- 
dependent substrate binding, (0ti) — ^x, is the mechanism 
behind this relation. 




Fig. 5: The distributions of (^tr) — ^btm (a) and (6'tr) — 0^ (b). 
36 molecules (108 states), c. The torque generation by the 
Fi-motor. The heat through the rotations, Q, is calculated as 
the potential change between successive transitions. When the 
transition is limited around Sx, A/i ~ Q. 

Discussions. — The 120° step triggered by the ATP 
binding is known to comprise 80° and 40° sub steps with a 
short pause after the 80° sub step, in which ATP hydroly- 
sis and phosphate release occur [TUlIIHlIin] . The duration 
of this pause is too short to be resolved by the relatively 
large probes used here. Therefore, the recovered potential 
is considered to be the effective one that the substate's 
potential is superimposed [20l[21]. In other words, with- 
out knowing the details of the system, we can recover the 
coarse-grained potentials to the extent of the model de- 
scription. It is the next challenge to use an essentially 
friction-free probe such as a colloidal gold particle with 
tens of nanometre diameter [TH1I21] for better estimation 
as well as the recovery of the potential profiles during the 
catalytic dwells. 

Recent developments in experimental techniques have 
enabled us to obtain single-molecule trajectories with high 
resolutions in both time and space. However, the tra- 
jectory analyses have thus far been limited to quantities 
such as the mean velocity and dwell time distribution; 
nevertheless, we expect that more valuable information 
is embedded in the trajectory. Despite the fact that our 
method requires only the single-molecule trajectory and 
the assumptions of the balance condition for the probe's 
motion in each chemical state, the Markovian proper- 
ties, and the number of states, it provides information 
about the state-specific potential profile, mechanochemi- 
cal coupling, and energetics. The method demonstrated 
that the mechanochemical picture of the torque gener- 
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ation is elucidated in terms of the basic concepts such 
as "power stroke" and "diffusion catch" and that the 
position-dependent substrate binding is the key feature of 
the mechanism behind the efficient free-energy conversion 
of the Fi-motor. The method is simple and will provide a 
promising methodology for understanding the properties 
of not only molecular motors but also systems in a broad 
range of fields. 

Materials and methods. — The experimental sys- 
tem was essentially the same as the previous studies [SJE] . 
Rotation of Fi-motor molecules (03/337 sub complex) de- 
rived from a thermophilic Bacillus PS3 was probed by at- 
taching strcptavidin-coated polystyrene particles (diame- 
ter = 0.3 /im, Seradyn) to the biotylated 7-shaft. Single- 
molecule observation was performed in the presence of 0.5 
MgATP, 2.5 MgADP, and 1 mM P; in a buffer 
containing 5 mM MOPS-K and 1 mM MgCl2 on a mi- 
croscope (Olympus) at 9,000 Hz (exposure time = 40 ^s) 
using a lOOx objective (NA = 1.49, Olympus), a Xenon 
lamp, and a high-speed camera (Easier) at a room tem- 
perature (25 ± 1 °C). The trajectory consists of approxi- 
mately 2,500,000 frames for each molecule typically. The 
mean rotational velocity and the diffusion coefficient were 
1.9 ± 0.9 Hz and 14.3 ± 3.6 rad^/s, respectively (mean ± 
SD, 36 molecules). The data including long pause pre- 
sumably due to the MgADP-inhibited state are excluded 
from the analysis. Bin width for calculating Ui was 1°. 
The value of the standard free energy change of an ATP 
hydrolysis (A^°) depends on literatures. We took the av- 
erage of calculated on the basis of three references 
[23H25] and calculated Afi 

A, = A,° + fcBTln^J^. (6) 

* * * 
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